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ABSTRACT 

We investigate the non-gaussian properties of cosmic-string-seeded linear density per- 
turbations with cold and hot dark matter backgrounds, using high-resolution numerical 
simulations. We compute the one-point probability density function of the resulting 
density field, its skewness, kurtosis, and genus curves for different smoothing scales. A 
semi-analytic model is then invoked to provide a physical interpretation of our results. 
We conclude that on scales smaller than 1.5(Qh 2 ) Mpc, perturbations seeded by cosmic 
strings are very non-gaussian. These scales may still be in a linear or mildly non-linear 
regime in an open or A- universe with V = flh < 0.2. 

Subject headings: cosmic strings — dark matter — galaxies: clusters: general — large- 
scale structure of universe 
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1. Introduction 

At present, the two main candidates for the ori- 
gin of cosmic structure are inflation and topological 
defects (for a review, see Vilenkin & Shellard 1994). 
Although both scenarios may produce a power spec- 
trum of density perturbations consistent with obser- 
vations, they have very different predictions regarding 
the statistical properties of the density field. While 
most inflationary models produce gaussian random- 
phase initial conditions, defect models produce non- 
gaussian perturbations particularly on small scales. 
New results from cosmic-string-seeded structure for- 
mation using high-resolution simulations (Avelino et 
al. 1997, 1998) were encouraging for models with 
r = Qh = 0.1-0.2 (see also Battye et al. 1997); both 
the mass fluctuation amplitude at 8/i _1 Mpc, as, and 
the power spectrum shape of cosmic-string-induced 
cold dark matter fluctuations, V(k), were consistent 
within uncertainties with observational data (Peacock 
& Dodds 1994; Viana & Liddle 1996). However, be- 
cause cosmic strings induce non-gaussian density per- 
turbations on small scales, the power spectrum alone 
is insufficient to describe all the statistical properties 
of such a density field. This is even more important in 
open or A-models because in those models the charac- 
teristic scales of the density field are shifted to larger 
scales relative to a fiat model with A = 0. 

In this Letter we investigate the non-gaussian prop- 
erties of the linear density field induced by cosmic 
strings using higher-order statistics such as the skew- 
ness and the kurtosis of a one-point probability den- 
sity function (PDF), as well as genus statistics. The 
non-gaussian properties we reveal provide a signifi- 
cant observational signature for cosmic string-seeded 
structure formation models on small length-scales. 
We note that previous analytic work has investi- 
gated the string-induced velocity field on scales above 
several h~ 1 Mpc, which was inferred to be gaussian 
(Vachaspati 1992; Moessner 1995), and that some of 
the features we study here were also observed in global 
topological defect models, notably for textures (Park, 
Spergel, & Turok 1991). Past work on genus statis- 
tics in the context of topological defects was made 
using toy models which incorporated some important 
features of the models in question (Brandenberger, 
Kaplan & Ramsey 1993; Albrecht & Robinson 1995; 
Avelino 1997). 

Our first step in the present analysis was to per- 
form high-resolution numerical simulations of cosmic 



string networks in an expanding universe (Allen & 
Shellard 1990) from which we subsequently computed 
the causally-sourced density field with either a cold or 
hot dark matter (CDM or HDM) background. The 
cosmic string simulations had a dynamic range ex- 
tending from before the radiation-matter transition 
at 0.4?7 e q through to deep into the matter era 8.4r/ cq , 
where r] cq is the conformal time at radiation-matter 
density equality. The structure formation simulation 
boxes contained 256 3 grid-points and their physical 
volume was in the range (4-100/i _1 Mpc) 3 . A much 
more detailed description of these methods is given 
by Avelino et al. (1997, 1998). 

2. Non-gaussian test statistics 

We first convolved the density field obtained above 
with a top-hat window at different smoothing scales 
R, and then calculated the resulting PDF p{v), where 
v is the number of standard deviations from the mean. 
The chosen smoothing scales were always at least 
three times larger than the grid spacing to guarantee 
a sufficient resolution for the smoothing, and at least 
twelve times smaller than the box size to reduce sam- 
ple variance effects. We then used the skewness and 
the kurtosis to test their gaussianity. The skewness is 
defined as M 3 /M^ 2 , where Mj is the i-th central mo- 
ment of a PDF. Its sign indicates the skewed direction 
against a normal distribution. The kurtosis is defined 
as M4/M2 — 3 and measures the size of the side tails 
against those of a normal distribution. Both statistics 
are zero for a gaussian distribution, but will depart 
from zero when the distribution is non-gaussian. We 
can therefore estimate a "non-gaussian scale" i? ng , be- 
low which both the skewness and the kurtosis depart 
significantly from zero. 

Another useful statistic sensitive to the non-gaus- 
sian features of a density field is the genus, which 
is employed to measure the topology of a continuous 
field. Conceptually, the genus of a given surface can 
be defined as g = # of holes — # of isolated regions + 
1. A more useful definition of genus is given by the 
Gauss-Bonnet theorem which relates the integrated 
Gaussian curvature (a local property) of a surface 
with the genus (a global property) of that surface: 

j>KdA = ^{l-g), (1) 

where dA is a differential two-dimensional surface 
patch, and K = \ja\a<2, is the Gaussian curvature 
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Fig. 1. — The one-point PDF p(y) of the density field 
in the linear string+CDM (left) and string+HDM 
(right) models for different R (in (f2/i 2 ) _1 Mpc). Here 
(74c : eric : Co. 5c : C4h : fih : co.5h * 1 : 5.2 : 8.9 : 
0.69 : 1.2 : 1.3. All simulation curves are averaged 
from 6 realizations, with standard deviation less then 
5% of the amplitudes throughout. 



of the patch with ai and et 2 the two principal radii 
of curvature. We used the numerical code developed 
by Avelino (1997) to calculate the genus of these iso- 
density surfaces, applying methods proposed by Gott, 
Mellot & Dickinson (1986). Here we employed the 
commonly used volume fraction parametrization of 
the genus curve, that is, we define a new variable vf 
so that 

f = - ± TJ^I e-" 2 /*du, (2) 



(27T) 



Vf 



where / is a previously calculated volume fraction of 
the whole simulation box within which the density 
is above a given threshold. Having done this, the 
genus curve becomes less-sensitive to the PDF which 
we studied separately. 

3. Results 

In Figure |l| we plot the one-point PDF p(y) of the 
density field in the linear string CDM and HDM mod- 
els for different smoothing scales R. The skewness 
and kurtosis of p{y) as functions of R were plotted 
in Figure 0. Although on large scales string pertur- 
bations are almost indistinguishable from gaussian 
random-phase perturbations, on scales smaller than 
1.5(Q/i 2 ) Mpc the non-gaussian character is remark- 
ably distinct, especially in a CDM background. 

In Figure [|we plot the genus curves g(ff) for R = 
0.4(fi/i 2 ) _1 Mpc in both string+CDM and string+HDM 
models, as well as those of gaussian random fields 
with correspondingly identical power spectra. These 
curves have been smoothed using three-point box- 
car smoothing (see for example Vogeley et al. 1994) 
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Fig. 2. — The skewness and kurtosis of p(y) induced 
by cosmic strings. 



Gaussian — Gaussian 

Cosmic strings 120 Cosmic strings 




Fig. 3. — The genus curves in the CDM (left) and 
HDM(right) models. Here R = 0.4(f2/i 2 ) -1 Mpc and 
the box size is (50/i _1 Mpc) 3 . Error bars are taken 
from 6 realizations for each curve. 



and were parametrized by the volume fraction (g). 
We can observe that on these small scales the genus 
curves of the cosmic-string-seeded density fields devi- 
ate in a significant way from gaussian random fields. 
As expected, the correlations among the phases of 
different Fourier modes imply a smaller number of 
independent regions in real space and therefore a 
smaller genus amplitude when compared with that 
of a gaussian random field. This is illustrated in Fig- 
ure ||, where two different isodensity thresholds for the 
CDM string model are directly compared with purely 
gaussian fluctuations for R = 0.4(fi/i 2 )" 1 Mpc. How- 
ever, we have also verified that when smoothed on 
larger scales R > 3(0/i 2 ) Mpc, cosmic-string and 
gaussian-random-phase genus curves are very close 
to each other in both amplitude and shape. Fi- 
nally, we note that the string+HDM model depar- 
ture from a gaussian distribution is less apparent than 
for string+CDM model throughout all the analysis 
above. 
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Fig. 4. — The isodensity surfaces with Vf = 
(A, B) and vt = 2 (C, D) for a gaussian field (A, 
C) and the string+CDM (B, D) model. Here R = 
0.4(f2ft 2 ) -1 Mpc and the box size is (25/i~ 1 Mpc) 3 . 



4. Discussion 

In order to interpret our results, we invoked the 
semi-analytic model described in Avelino et al. (1997, 
1998; see also Albrecht & Stebbins 1992). This model 
was shown to be very accurate in reproducing the 
simulation power spectrum in both matter-era and 
radiation-era scaling regimes, as well as during the 
radiation-matter-era transition. The power spectrum 
of string-induced density perturbations can be written 
as: 

pit _ 

V(k) = 16^ 2 GV / \G{k; m, ??)l 2 ^, v)dv > (3) 

where /i is the string energy density per unit length, 
Q(k;rio,r)) the Fourier transform of the appropriate 
Green's functions, and T{k,rf) the structure func- 
tion. At a given time, J-(k,rj) has a turn-over scale 
at » 20/t; which reflects the correlation length of 
cosmic strings £ ~ At a particular time, the 

perturbations induced on scales larger than the string 
correlation length are generated by many string ele- 
ments and, therefore, they are nearly gaussian accord- 
ing to the central limit theorem. On the other hand, 
perturbations induced on smaller scales are very non- 
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Fig. 5. — The estimated gaussian and non-gaussian 
contributions as a comparison with the total power 
spectrum in the CDM (left) and HDM (right) models. 

gaussian because they can be either very large within 
the string-induced wakes or else very small in re- 
gions outside these wakes. In consequence, we can 
roughly divide the power spectrum of cosmic-string- 
seeded density perturbations into two parts: a nearly 
gaussian part generated when the string correlation 
length was smaller than the scale under considera- 
tion, and a strongly skewed non-gaussian part gen- 
erated when the string correlation length was larger. 
We shall abuse strict definitions and call these sim- 
ply the 'gaussian' and 'non-gaussian' contributions re- 
spectively. In terms of the structure function in (||) , 
we can make the split Fgik, rf) = J-(k,rf) for k < k^ 
(JF g = for k > k{) and, similarly, T n g{k,rj) = 
!F{k,rj) for k > k%. J- g corresponds to the left- 
hand side of the structure function from its peak to 
the causal compensation cut-off where J- g oc fc 4 for 
k « kf. J- ng can be largely identified with the im- 
print of string wakes and behaves as J- ng oc k~ 2 for 
k > fcg. 

Integrating (^) with this gaussian/non-gaussian 
split, we can compare the relative contributions to 
the total power spectrum in CDM and HDM mod- 
els, as illustrated in Figure H. This plot demon- 
strates that these two contributions become com- 
parable at k ~ 2/iMpc -1 for the CDM model and 
k ~ 0.5/iMpc" 1 for the HDM model, taking a Hubble 
parameter Hq — 70km s _1 Mpc _1 . We further calcu- 
lated the variance of the density fields from these non- 
gaussian and gaussian contributions when smoothed 
with a top-hat window of radius R: 

poo 

a 2 R = An \w(kR)\ 2 P(k)k 2 dk 7 (4) 
Jo 

where w(x) = 3ji(a;)/x is the Fourier transform of a 
top-hat window and ji(x) the spherical Bessel func- 
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Fig. 6. — The ratio crR(ng) /prm against R for both 
CDM (left) and HDM (right) models. The solid lines 
have a full dynamic range going from r\ — to today, 
while the dashed lines have the same dynamic range 
of our simulations. 

tion of order one. Figure || plots the ratio between 
the standard deviations of the non-gaussian and gaus- 
sian density perturbations, a _R( ng ) / 'a rm , against R 
for both CDM and HDM models. The criterion 
a R(ng) / & R(g) > lj gives an estimate for the non- 
gaussian scale i? ng - In Figure [| we also consider 
the effect of the limited dynamic range of our sim- 
ulations, by comparing with the full semi-analytic 
integral (^). The simulations slightly enhance non- 
gaussianity essentially because there is more 'gaus- 
sian' power missing from early times rj < 0.4?7 oq . 
Nevertheless, we have checked that the power lost 

(c(full) ~ o-( 8 im))/o-(full) on scales R - l(^ 2 ) -1 Mpc 
is always < 10% for the non-gaussian contribution, 
< 35% for the gaussian contribution, and < 25% for 
the total power. This confirms the reliability of our 
non-gaussian analysis based on these simulations. 

The consistency between Figure || and Figure || 
should be evident, notably for the skewness. The 
physical reason for the much stronger non-gaussian 
features on small scales can now be interpreted as the 
following. Before the radiation-matter equality rj eq , 
the growth of density perturbations was suppressed. 
This means that even on small scales the contribution 
to the power spectrum from the matter era is very im- 
portant. Adding the fact that by the time of r] eq the 
string correlation length had already grown to several 
(fi/i ) _1 Mpc, the perturbations generated during the 
matter era on scales smaller than this will add an im- 
portant non-gaussian contribution to the total power 
spectrum. This also explains the peak in the PDF 
on small scales (see Figure ^). The limited dynamic 
range of our simulations tends to produce a sharp 
peak on the PDF for small smoothing scales because 
of the presence of voids, that is, regions left largely 



unperturbed by the transition era strings. However, 
even for the full dynamic range a similar peak, though 
slightly smoother, should still emerge. This is the ori- 
gin of the non-gaussianity on smaller scales, the exis- 
tence of which is not affected by the limited dynamic 
range as we demonstrate in Figure | The reason 
why the HDM model is more gaussian than the CDM 
model is that the neutrino free-streaming length pre- 
vents small-scale perturbations from emerging until 
late times. This significantly reduces the HDM power 
at large k, precisely where the non-gaussian contribu- 
tion is dominant in the CDM model (see Figure |j). 
This also accounts for the smaller genus amplitude for 
the HDM model (see Figure §). 

We also note that our results were obtained in the 
linear regime. When the density field evolves into 
the non-linear regime at later times, the skewness in- 
creases as a direct result of the constraint that the 
density contrast 5 > — 1. However, when compared 
with observations, the skewness of our linear cosmic- 
string-seeded result for any reasonable value of T is 
well below that of the QDOT survey (Saunders et 
al. 1991). Our results are also consistent with those 
of Canavezes et al. (1997) as we find that cosmic- 
string and gaussian-random-phase genus curves are 
very similar in the scale range where no strong phase 
correlations were found by the genus analysis of the 
PSCz survey. For gaussian models, it is well known 
that in the non-linear regime the skewness is equal to 
its initial value plus a contribution due to the non- 
linear collapse. This contribution scales proportion- 
ally to the rms density fluctuation a (Fry & Scherrer 
1994). Similarly, the non- linear contribution to the 
kurtosis also has a simple scaling dependence on er, 
depending on whether or not the initial density field 
is gaussian (Chodorowski & Bouchet 1995). This in- 
dicates that non-linear clustering preserves the non- 
gaussianity inherent in the linear evolution, but fur- 
ther investigation using iV-body simulations will be 
required to establish if the non-gausian signature of 
cosmic string models is potentially measurable. 

5. Conclusion 

We conclude that on length scales smaller than 
l.5(tth 2 ) 1 Mpc perturbations seeded by cosmic strings 
are very non-gaussian, especially in the context of a 
CDM model. In an open or A-universe with T — Vlh ~ 
0.15, this scale will be shifted to 10/i~ 1 Mpc, which 
may still be in the linear or mildly non-linear regime, 
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thus potentially providing a strong empirical test for 
cosmic string models. It has been suggested that such 
non-gaussianity may imply that it is difficult in string 
models to deduce the parameter (3 = S1q' 6 /6 from ob- 
servations of density and velocity fields (van de Bruck, 
1997). However, our results indicate otherwise on 
large scales because we find that string perturbations 
are very similar to gaussian-random phase fluctua- 
tions, when smoothed on sufficiently large scales. 
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